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Resistive magnetohydrodynamic simulations of X-line retreat during 
magnetic reconnection 
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To investigate the impact of current sheet motion on the reconnection process, we perform resistive magne- 
tohydrodynamic (MHD) simulations of two closely located reconnection sites which move apart from each 
other as reconnection develops. This simulation develops less quickly than an otherwise equivalent single 
perturbation simulation but eventually exhibits a higher reconnection rate. The unobstructed outflow jets 
are faster and longer than the outflow jets directed towards the magnetic island that forms between the two 
current sheets. The X-line and flow stagnation point are located near the trailing end of each current sheet 
very close to the obstructed exit. The speed of X-line retreat ranges from ~0. 02-0. 06 while the speed of 
stagnation point retreat ranges from ~0. 03-0. 07, in units of the initial upstream Alfven velocity. Early in 
time, the flow stagnation point is located closer to the center of the current sheet than the X-line, but later 
on the relative positions of these two points switch. Consequently, late in time there is significant plasma 
flow across the X-line in the opposite direction of X-line retreat. Throughout the simulation, the velocity at 
the X-line does not equal the velocity of the X-line. Motivated by these results, an expression for the rate of 
X-line retreat is derived in terms of local parameters at the X-point. This expression shows that X-line retreat 
is due to both advection by the bulk plasma flow and diffusion of the normal component of the magnetic field. 

PACS numbers: 52.35.Vd,52.65.-y 

Keywords: Magnetic reconnection, plasma simulation, plasma magnetohydrodynamics 



I. INTRODUCTION 

Most simulations and theories of magnetic reconnec- 
tion operate under the assumption that the current sheet 
is roughly stationary with respect to the ambient plasma. 
However, there are many situations in nature and the 
laboratory where current sheet motion is important. Re- 
cently, a model was developed to describe steady mag- 
netic reconnection with asymmetric outflow in a high as- 
pect ratio current sheet. 1 While the assumption of time- 
independence was needed to make analytic progress, such 
an assumption precludes time-dependent effects such as 
current sheet motion. This paper addresses this issue by 
presenting resistive magnetohydrodynamic (MHD) sim- 
ulations of two X-lines which start in close proximity to 
each other and move apart as reconnection develops. 

The near-Earth neutral line model 2 predicts that the 
magnetotail X-line retreats in the tailward direction dur- 
ing the recovery phase of magnetospheric substorms. 
Such behavior is commonly observed during in situ mea- 
surements in the Earth's magnetotail^ In a recent sta- 
tistical study of diffusion region crossings by Cluster, 
tailward moving X-lines were observed ~4 times more 
frequently than earthward moving X-lines^ Despite the 
common occurrence of current sheet motion, it is stan- 
dard practice to compare in situ measurements of diffu- 
sion region crossings to particle-in-cell (PIC) simulations 
of roughly stationary reconnection layers. Global simu- 
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lations of the magnetotail do allow the diffusion region 
to move, and tailward motion of the predominant X-line 
is commonly observed A"— Previous research has also con- 
sidered the effects of current sheet motion on reconnec- 
tion slow mode shock structured- 
Flux rope models of coronal mass ejections (CMEs) 
typically predict the formation of a current sheet be- 
tween the flare site and the ejected plasmoid.— ~— Fea- 
tures identified as current sheets have been observed for 
an increasing number of events (e.g., Refs. Fl5l - fl8h . Dur- 
ing these events both the lower and upper boundaries of 
the current sheet are thought to rise with time i 19 i 2Q Of 
particular interest are the locations of the predominant 
X-line and flow stagnation point. A recent analysis of the 
current sheet behind a slow CME on 2008 April 9 showed 
upflowing features above and downflowing features be- 
low a height of ~0.25 solar radii. 18 In addition, a po- 
tential field model of the pre- CME active region showed 
a pre-existing X-line at about that height. Because the 
CME current sheet associated with this event extended 
beyond several solar radii, this evidence suggests that the 
predominant X-line and flow stagnation point are both 
located near the base of the current sheet. 

Current sheet motion and asymmetric outflow recon- 
nection occur in laboratory plasma devices involving the 
merging of spheromaks and toroidal plasma configura- 
tions where the outflow is aligned with the radial di- 
rection, even though the range of motion is limited by 
the boundary conditions of the experiment. Counter- 
helicity spheromak merging results from the Magnetic 
Reconnection Experiment (MRX) 21 show the X-line be- 
ing pulled towards one end of the current sheet because of 
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the Hall effect, resulting in asymmetric outflow j 22 i 23 Dur- 
ing spheromak merging experiments at TS-3/4, 24 cur- 
rent sheet ejection often leads to faster reconnect ion and 
the X-line being located towards one end of the current 
sheets 

Current sheet motion also occurs when both a guide 
field and a density gradient across the current layer are 
present^"— Relevant configurations include the dayside 
magnetopause and tokamaks. In these situations the 
plasma pressure gradient in the inflow direction leads 
to diamagnetic drifting of the reconnect ion layer. When 
the diamagnetic drift velocity becomes comparable to the 
Alfven speed, reconnection is suppressed. Such effects 
can contribute to mode rotation in tokamaks. 

Recent PIC simulations by Oka et al 29 displayed X- 
line retreat due to the presence of an obstructing wall in 
one of the two downstream regions. The retreat speed 
was found to be ~0.1 of the upstream Alfven velocity 
and comparable to the reconnection inflow velocity. In 
qualitative agreement with the scaling model developed 
in Ref. 1, the reconnection rate was not significantly af- 
fected by current sheet motion. However, the structure of 
the current sheet was modified from the symmetric case. 
In particular, the ion flow stagnation point was closer to 
the wall than the X-line, and the outflow jets away from 
the obstruction were longer and faster than the outflow 
jets directed towards the wall. 

In all of these cases, current sheet motion is an im- 
portant consideration because it modifies the structure 
of the reconnection layer and affects the transport of en- 
ergy, mass, and momentum. The simulations in this pa- 
per are used to help understand the impact current sheet 
motion has on the structure and dynamics of a reconnec- 
tion layer and to determine what sets the rate of X-line 
retreat. 

This paper is organized as follows. Section Hfl describes 
the numerical method used for the simulations reported 
in this paper. Section HIJ describes the initial perturbed 
equilibrium and boundary conditions. Section [IV] pro- 
vides a thorough discussion of the simulation results, in- 
cluding comparisons to a symmetric single perturbation 
simulation and details of the internal structure of the cur- 
rent sheet. Section IVl contains a derivation of the rate of 
X-line retreat in terms of local parameters evaluated at 
the X-point. Section [VTl contains a summary and conclu- 
sions. 



II. NUMERICAL METHOD 

The NIMROD code^ 2 . (Non-Ideal Magnetohydrody- 
namics with Rotation, Open Discussion) is well suited 
for the study of resistive MHD and two-fluid magnetic 
reconnection in a variety of configurations j 23 i 32 i 33 NIM- 
ROD uses a finite element formulation for the poloidal 
plane and, for three dimensional simulations, a finite 
Fourier series in the out-of-plane direction. The equa- 
tions evolved by NIMROD for the two-dimensional sim- 



ulations reported in this paper are given in dimensionless 
form by 

% + V • (pV) = V • DVp, (1) 

— = -Vx(^J-VxB), (2) 
J = VxB, (3) 

p + V • VV) = JxB-Vp + V- pvW, (4) 

^l(f + v.vr) = -fv.v-v. q + Q, (5) 

where B is the magnetic field, V is the bulk plasma veloc- 
ity, J is current density, p is density, p is the plasma pres- 
sure, r] is resistivity, v is the kinematic viscosity, T is tem- 
perature, q = — pxVT represents isotropic thermal con- 
duction, x is the thermal diffusivity, Q includes resistive 
and viscous heating, 7 = 5/3 is the ratio of specific heats, 
and D is an artificial number density diffusivity. Simu- 
lation quantities are normalized to the following respec- 
tive values: B , p , L , t > Vao = L /t = 5 /v7^oPo> 
Po PoV^o B /fi poTo/rrii, J B /poLo, and 
W^o ^0 = Xo = A) = Ll/t . The divergence con- 
straint is not exactly met, so divergence cleaning is used 
to minimize the development of divergence error. 30 

NIMROD represents solution fields as the sum of 
steady-state and time- varying components. For most ap- 
plications the use of an ideal MHD equilibrium for the 
steady-state component does not lead to significant error 
when a small finite resistivity is used. However, in this 
work it is necessary for the steady-state component to 
be free of current for resistive diffusion to be represented 
accurately. 



III. PROBLEM SETUP 

The initial conditions for the two-dimensional simu- 
lations reported in this paper are those of a perturbed 
Harris sheet equilibrium. The Harris sheet equilibrium is 
given by 



B x (z) = Bo tanh 
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where Bo is the asymptotic magnetic field strength, £0 
is the Harris sheet thickness, and /3q is the ratio of the 
asymptotic plasma pressure to the asymptotic magnetic 
pressure. The initial temperature is constant. Here, x 
is the outflow direction, y is the out-of-plane direction, 
and z is the inflow direction. The perturbed component 
of the magnetic field is of the form 



B p 0, z) = V x (V^y) 



(9) 
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with the flux function tp p for double perturbation simu- 
lations given by 



ip P (x,z) = - 




where B\ is the amplitude of the perturbation, A is the 
half-separation between the two components of the per- 
turbation, and h is the length scale of the perturbations. 
This configuration is reminiscent of the final state of the 
multiple X-line simulations presented in Ref. 34. The 
form of this perturbation does not depend on the size of 
the computational domain. 

The computational domain is defined by — x max < x < 
x max and -Zmax < z < z max . The domain is periodic 
in the x direction and has no-slip perfectly conducting 
boundaries at z = ±z max . The domain has m x and m z 
finite elements in the x and z directions, respectively. 
Mesh packing is used in both the x and z directions to 
ensure that the reconnection layers are sufficiently re- 
solved. The physical size of the domain is chosen to 
be large enough that boundary conditions do not signifi- 
cantly affect the dynamical behavior over the time scales 
considered in the analysis. 

The simulation parameters are as follows. The Har- 
ris sheet equilibrium is given by Bq = 1, Sq — 0.1, and 
Po = 1. The initial perturbation is given by B\ = 0.05, 
A = 1, and h = 0.5. The length scale is normalized to the 
half-distance between the initial perturbations. The dif- 
fusivities are given by r] = v = x = 10 -3 and D = 10 -5 . 
All velocities given in this paper are in the simulation 
reference frame. A variable time step is used such that 
the Courant number is no greater than 1/6. The com- 
putational domain is given by x max = 30, z max = 12, 
m x = 124, and m z — 22 with sixth order finite element 
basis functions. The domain is chosen to be large enough 
that boundary conditions do not significantly affect the 
relevant behavior for the time scales considered in this 
paper. The simulations are ended when the outflow jets 
reach the periodic boundary. Because the simulation is 
symmetric about x = 0, only x > is considered in the 
following sections. 



IV. SIMULATION RESULTS 

This section contains a detailed analysis of a double 
perturbation simulation which is used to understand the 
structure and development of a moving current sheet in 
the resistive MHD limit. This double perturbation simu- 
lation is compared to an otherwise equivalent single per- 
turbation simulation which exhibits non-retreating, sym- 
metric reconnection. 



A. The global structure of the reconnection region 

Key features of the time evolution of the double pertur- 
bation simulation are shown in Fig. [TJ The most appar- 
ent feature of the simulations is that the current sheets 
retreat from x = while the reconnection process is de- 
veloping. The current sheet has a single wedge shape 
which is most apparent late in time. Figure [2fa) shows 
that the current sheet is thinnest near the magnetic field 
null, which is very close to the obstructed exit from the 
current sheet, and thickest near the obstructed exit. 

As seen in Figs. [TJd)-(f), the peak plasma pressure in 
the magnetic island around x = is greater than the peak 
plasma pressure in the downstream region away from 
x = 0. The plasma pressure gradient associated with 
the obstructing magnetic island is stronger and more lo- 
calized than the plasma pressure gradient associated with 
the unobstructed outflow region. Early in time (t < 20), 
the X-point is located near the minimum of plasma pres- 
sure along z = 0. Later (20 < t < 28), the plasma 
pressure at the flow stagnation point is less than at the 
X-point. A local maximum of plasma pressure appears 
at t « 27 in the central region of the current sheet which 
persists through the remainder of the simulation. This 
pressure maximum is located slightly closer to the un- 
obstructed exit of the current sheet than the obstructed 
exit. 

Contours of the outflow component of velocity pre- 
sented in Figs.[]Jg)-(i) show that the outflow jet directed 
away from the obstruction is faster and longer than the 
outflow jet directed towards the obstruction. This behav- 
ior is consistent with previous simulations of asymmet- 
ric outflow reconnection£ i 2Q i 29 i 34 ~ — The kinetic energy 
and enthalpy fluxes away from the obstruction are much 
greater than the those towards the obstruction. 

The inflow component of velocity is shown in Figs. 
[Uj)-(i). In the upstream regions outside the current 
sheet, the inflow component of velocity is roughly uni- 
form except near the outflow region furthest from x = 0. 
The inflow velocity late in time is ~0. 01 5-0. 02. Because 
the large outflow blob from the unobstructed exit of the 
current sheet seen in Fig. [1] is propagating away from 
x = 0, there is a quadrupole-like pattern for V z as up- 
stream plasma ahead of the blob moves out of the way 
and upstream plasma behind the blob is pulled back to- 
wards z = 0. 

Figure [3] shows the time evolution of key quantities 
from the double perturbation simulation with compar- 
isons to the otherwise equivalent single perturbation sim- 
ulation used as a control for this study. Because a mag- 
netic island forms at t « 38 in the single perturbation 
simulation, the comparison is halted at that time. Fig- 
ure [3ja) shows the reconnection electric field strength 
for both simulations. The reconnection rate peaks more 
quickly in the single perturbation simulation than the 
double perturbation simulation. The development of re- 
connection in the double perturbation simulation is slow 
because the magnetic field null is located near the global 
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FIG. 1. (Color online) Contour plots showing: (a)-(c) the out-of-plane current density, J y , (d)-(f) the plasma pressure, p, 
(g)-(i) the outflow component of velocity, V x , and (j)— (1) the inflow component of velocity, V z , at times 22, 44, and 66 for one 
of the two reconnection regions during the simulation of X-line retreat discussed in Sec. lIVI The locations of the magnetic field 
null, flow stagnation point, pressure minimum along z = 0, and pressure maximum within the current sheet are denoted by 'n,' 
's,' 'pmin,' and 'pmax.cs,' respectively. 



pressure minimum along z = for t < 20 [see Fig. Efd)] 
so that pressure gradients oppose outflow instead of fa- 
cilitating it. However, the amplitude of the reconnection 
rate is eventually greater in the double perturbation case 
because the current sheet can increase in length in only 



one direction. As seen in Fig. [3fb), the single perturba- 
tion current sheet is a factor of ~2-3 longer at any given 
time. The slight disturbances around t > 58 in Figs. (Ha) 
and [31(b) may be due to the periodic boundary condition. 
While the reconnection rate is enhanced slightly in the 
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FIG. 2. (Color online) Shown above at t = 66 are (a) the 
current sheet thickness 6 as a function of x, calculated as 
the half width at quarter maximum of J y , and (b) the plasma 
pressure along z = 0. In (a), the magnetic field null is denoted 
by 'n' and the flow stagnation point is denoted by 's,' with 
x n = 3.47 and x s — 4.05. In (b), the global pressure minimum 
along z — and the pressure maximum within the current 
sheet are denoted by min ' and l v max,CS 5 ' respectively, with 
x p ,min = 2.98 and x p , ma x,cs = 9.90. 



double perturbation case, the difference is less than a 
factor of two so we conclude that the reconnection rate 
is only modestly affected by asymmetry or current sheet 
motion (see also Refs. U andl29h. 

The peak outflow velocities in the positive and negative 
directions are shown in Fig. [3jc). Again, it is apparent 
that the double perturbation simulation takes longer to 
develop than the single perturbation simulation. The 
unobstructed outflow jet is typically ~2-3 times faster 
than the obstructed outflow jet. 



B. The internal structure of the current sheet 

The positions of the magnetic field null, flow stagna- 
tion point, and pressure minimum along z = (x ni x s , 
and Xp 5 min, respectively) are presented in Fig. [3f d). The 
positions of the X-point and flow stagnation point are 
separated by a short distance. Both points are located 
very close to the obstructed exit of the current sheet. 
Early in time, the flow stagnation point is closer to the 
obstruction than the magnetic field null (e.g., x s < x n ). 
This corresponds to the plasma flow at the X-line moving 
in the same direction as X-line retreat. Such a scenario 
is what one would expect if the frozen-in condition is ap- 
proximately met [see Fig. Ufa)]. However, around t « 28 
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FIG. 3. (Color online) Simulation results showing as a func- 
tion of time (a) the reconnection electric field strength mea- 
sured at the X-point, (b) the full length at quarter maximum 
of J y along z = 0, (c) the peak outflow velocities, (d) the po- 
sitions of the magnetic field null, flow stagnation point, and 
pressure minimum along z = (x n , x s , and x p , m in, respec- 
tively), and (e) the rates of change in position of the magnetic 
field null and flow stagnation point (dx n /dt and dx s /dt, re- 
spectively) and the velocity at the magnetic field null, V x (x n ). 



the relative positions of the magnetic field null and flow 
stagnation point switch (e.g., x s > x n ). Hence at late 
times during this simulation the plasma flow at the X- 
line is in the opposite direction of X-line retreat [see Fig. 
2Kb)]. The peak resistive electric field in the double per- 
turbation simulation occurs slightly after the time when 
(e.g., near the time when tension forces entirely 
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(a) Early in time 




(b) Late in time 




FIG. 4. (Color online) A schematic showing the rate of X- 
line retreat, dx n /d£, and the plasma flow velocity at the X- 
line, V x (x n ). Early in time, both velocities are in the same 
direction. Late in time, there is significant plasma flow across 
the X-line in the opposite direction of X-line retreat. 



work towards accelerating outflow rather than partially 
working against it). 

The relative positions of the X-line and flow stagna- 
tion point switch because during near steady conditions, 
the flow stagnation point will in general be located near 
where the tension and pressure gradient forces cancel 
(e.g., Ref. 1). This implies that the total pressure at 
the flow stagnation point should be greater than the to- 
tal pressure at the X-point when time-dependent effects 
are not important. Early in time (t < 20), the X-point 
is located very near the global pressure minimum along 
z = [see Fig. [3f d)] . During this phase, p(x s ) > p(x n ) 
and the pressure gradient and tension forces at the flow 
stagnation point are oppositely directed. As the simu- 
lation progresses, the X-line and flow stagnation point 
retreat more quickly than the position of this pressure 
minimum. For 20 < t < 28, p(x s ) < p(x n ) and the ten- 
sion and pressure gradient forces at the flow stagnation 
point are pointed in the same direction and thus cannot 
cancel each other out. The flow stagnation point retreats 
quickly until x n < x s and p(x s ) > p(x n ) once again and 
for the remainder of the simulation [see Fig. EJb)]. 

Figure [3^e) compares the rate of change in position of 
the magnetic field null dx n /dt, the rate of change in posi- 
tion of the flow stagnation point dx s /dt, and the plasma 
flow velocity across the magnetic field null V x (x n ). As 
discussed in Ref. |37| and Sec. [VJ any difference between 
V x (x n ) and dx n /dt must be due to resistive diffusion. 
The X-line is observed to retreat at a variable veloc- 
ity ranging between 0.02 and 0.06. The flow stagnation 
point retreats more quickly than the magnetic field null 
for most of the simulation with velocities ranging from 
0.03 to 0.07. For much of the first third of the simulation, 
dx s /dt ~ dx n /dt ~ V x (x n ). When x s « x n , the rate of 
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FIG. 5. (Color online) The normal component of the magnetic 
field as a function of z for three locations surrounding the 
magnetic field null at t = 66. At x = x n , B z < in the 
vicinity of the current sheet except at z = 0. 
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FIG. 6. (Color online) Magnetic flux contours near the mag- 
netic field null at t — 66. The dotted line represents X — Xn • 



X-line retreat drops by a factor of ~3 before slowly in- 
creasing and stabilizing. Following a delay after the drop 
in dx n /dt : dx s /dt declines also before stabilizing. How- 
ever, the most striking feature of Fig. [3fe) is that late in 
time there is significant plasma flow across the X-line in 
the opposite direction of X-line retreat: V x (x n ) ~ —0.16 
while dx n /dt - 0.04. 

Figure [5] shows B z (z) for three locations surrounding 
the magnetic field null. Along x = x n , B z < except 
at z = 0. Additionally, d 2 B z /dt 2 < in the vicinity 
of the X-point. This magnetic field structure near the 
X-point can be seen clearly in the magnetic flux contour 
plot presented in Fig. [6j where magnetic field lines are 
pinched inward towards the right of the X-point. As we 
shall see in Sec. El these features facilitate X-line retreat 
by diffusion of the normal component of the magnetic 
field in the inflow direction. 

Plasma flow at the X-line in the direction opposite to 
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Inside the current sheet (3 



< <r. < 



14), the electric 



field is approximately uniform but increases slightly with 
distance from x = 0. The resistive electric field peaks in 
between the magnetic field null and flow stagnation point. 
In the obstructed outflow region, the electric field drops 
to a fraction of the value from inside the current sheet. 
From Faraday's law, the positive slope of E y indicates 
that B z in the obstructing magnetic island is increasing 
with time. In the unobstructed outflow region, E y peaks 
sharply. This signature is indicative of the unobstructed 
downstream plasmoid being advected away from x = 0. 

Figure E^b) shows force density contributions from the 
momentum equation (Eq. 2J along z = at t = 66. Con- 
tributions include the magnetic tension force B z dB x /dz, 
the magnetic pressure gradient force —d(B^/2)/dx, and 
the plasma pressure gradient force —dp/dx. Viscous 
forces are small and not shown explicitly, but are included 
in the calculation for total force density. The outflow to- 
wards the unobstructed exit is accelerated primarily by 
magnetic tension, whereas the outflow towards x = is 
accelerated primarily by the pressure gradient force. In- 
tuitively, this is because the X-line is located very near 
one end of the current sheet so that the tension force 
directed towards that end is small and the tension force 
directed away from that end is large (e.g., Ref. 23). The 
force density at the flow stagnation point is negative un- 
til t « 39 when it becomes positive for the remainder of 
the simulation. 



FIG. 7. (Color online) Simulation results showing (a) compo- 
nents of the electric field and (b) force density contributions 
along z = at t = 66. The magnetic field null is denoted by 
'n' and the flow stagnation point is denoted by 's.' 



X-line retreat appears to be a persistent feature for the 
assumed configuration. This effect occurs during analo- 
gous simulations with initial uniform guide fields up to 
at least B y o = 4, a line-tied boundary at x = 0, (3q 
at least between 0.25 and 4, different diffusivities, and 
different initial separations at least for A < 8. During 
some of these simulations, an additional X-line appears 
near the leading end of the current sheet. For larger A, 
the developing current sheets may become unstable to 
the plasmoid instability^ - — before the retreat process 
begins. Plasma flow across the X-line in the direction 
opposite to X-line retreat does occasionally occur during 
analogous resistive MHD simulations of multiple compet- 
ing reconnect ion sites, but less frequently because current 
sheet motion is inhibited. 44 



C. Components of the electric field and momentum 
balance 

Components of the out-of-plane electric field along z = 
are shown in Fig. [T^a) for t = 66. This is late in the 
simulation when reconnection is well-developed and x n < 



V. THE RATE OF X-LINE RETREAT 

This section contains a derivation of the rate of X-line 
retreat in the two-dimensional case. The geometry is as- 
sumed to be the same as that of the simulations reported 
in this paper. In particular, symmetry is assumed about 
z = so that the motion of the X-line is purely in the x 
direction. No assumptions are made about the presence 
or absence of a guide field, except that d/dy —> 0. The 
analysis presented in this section can also be applied to 
the asymmetric inflow case by switching the inflow and 
outflow directions (see also Ref. 1451) . 

The derivation of the rate of retreat of the X-line begins 
with Faraday's law, 



^ = - VxE - 



(ii) 



Evaluating the z-component of this expression and using 
the assumption of symmetry in the out-of-plane direction 
gives 



dB z 

at 



dx 



(12) 



Next, the convective derivative of B z at the X-point taken 
using the velocity of X-line retreat dx n /dt is given by 



dB z 



dt 



dx n dB z 



dt dx 



0. 



(13) 
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FIG. 8. (Color online) The mechanism for X-line retreat due 
to diffusion of the normal component of the magnetic field in 
the inflow direction. 



The right hand side of Eq. [13] is zero because the normal 
component of the magnetic field at the X-point does not 
change from zero, by definition. The rate of X-line retreat 
then follows from Eqs. [T2l and [T3l and is given by 



dx n 
dt 



dEy/dx 



dB z /dx 



(14) 



Intuitively, this means that the X-line moves in the di- 
rection of increasing reconnection electric field strength. 
This expression shows that it is not self-consistent to as- 
sume that the X-line is moving while the out-of-plane 
electric field is constant in space (cf. Ref. M). 

This derivation has thus far not specified an expression 
for the out-of-plane electric field. By using the resistive 
MHD Ohm's law E + VxB = 77J, Eq. [H yields 



because d 2 B z /dz 2 < 0, the normal component of the 
magnetic field becomes more negative in the immediate 
vicinity of the X-point. This diffusion of B z in the z di- 
rection causes the magnitude of B z to become stronger 
immediately to the left of the X-point and weaker im- 
mediately to the right of the X-point in Fig. [51 As a 
consequence, the X-point retreats to the right. The nec- 
essary features for this diffusion mechanism are apparent 
in Figs. [5] and [6l The contribution to X-line motion from 
diffusion of the normal component of the magnetic field 
in the outflow direction is not shown in Fig. [8] 

The difference between the bulk plasma flow across 
the X-line, V x (x n ), and the rate in change of position of 
the X-line, dx n /dt, is subtle but important. In the ideal 
MHD limit, the magnetic field will be purely frozen in to 
the plasma so that these two quantities will be identical 
[e.g., V x (x n ) = dx n /dt\. In this limit the X-line is purely 
advected by the bulk plasma flow. Therefore, any differ- 
ence between V x (x n ) and dx n /dt must be due to resistive 
diffusion of the magnetic field (see also Ref. I37I ). 

The relation presented in Eq.[T5land Fig. [8] shows that 
the rate of X-line retreat depends fundamentally on local 
parameters near the X-point. Hence global models at- 
tempting to describe the rate of X-line retreat must take 
into account the coupling between large scales and the 
local structure of the magnetic field near the X-point. A 
starting point for such models may be the analysis pre- 
sented in Ref. 0. 



VI. DISCUSSION AND CONCLUSIONS 



dXn 
dt 



V x (x n ) - T] 



d 2 B z 
dx 2 



2 B Z 



dz 2 



8B Z 
dx 



(15) 



For the assumed geometry, this is an exact result for resis- 
tive MHD. This analysis can be extended to include ad- 
ditional terms in the generalized Ohm's law. The inverse 
dependence on dB z /dx\ x in the resistive term comes 
from the geometric properties of Eq. [13j it is easier to 
change the position of a root of a function by a vertical 
shift if the slope of the function near the root is shallow. 
Because d 2 B z /dx 2 \ x will usually be small if dB z /dx\ Xn 
is small, it is likely that in resistive MHD the term 
d 2 B z /dz 2 \ representing diffusion of the normal com- 
ponent of the magnetic field in the inflow direction will be 
more important than the term d 2 B z / dx 2 1 representing 
diffusion of the normal component of the magnetic field 
along the outflow direction. Indeed, the simulation does 
show that late in time d 2 B z /dx 2 \ <C d 2 B z /dz 2 \ . 
When the X-line is retreating during collisionless recon- 
nection, it is possible that the term d 2 B z /dx 2 will be- 
come more important. During three-dimensional recon- 
nection, the term d 2 B z /dy 2 may also play a role in facil- 
itating X-line retreat. 

The mechanism for diffusion of the X-line is presented 
in Fig. [51 Along x = x ni B z < except at z = 0. Then 



This paper presents a resistive MHD simulation of two 
competing reconnection sites which move apart from each 
other as they develop. This investigation provides insight 
into the impact of current sheet motion on the reconnec- 
tion process as well as on what happens when outflow 
from one exit of the current sheet is blocked. 

The asymmetry in the reconnection process is most 
apparent in the outflow velocity profile. When the re- 
connection process is well-developed, the unobstructed 
outflow jet is ~2-3 times faster than the obstructed out- 
flow jet. As a consequence, most of the mass, energy, and 
momentum flux associated with the outflow is directed 
away from the obstructing magnetic island that forms 
between the two current sheets. Late in time, the recon- 
nection rate is slightly higher in the double perturbation 
simulation than in an otherwise equivalent single pertur- 
bation simulation because the obstruction prevents the 
length of each current sheet from growing in one direc- 
tion. 

The X-point and flow stagnation point are both located 
very near the obstructed exit of the current sheet. This 
gives the reconnection layer a characteristic single wedge 
shape which is especially apparent late in time. During 
recent simulations of the plasmoid instability^ - — single 
wedge shaped small-scale reconnection sites are routinely 
observed with the X-point located near the thinnest part 
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of the current sheet (see the online movie associated with 
Fig. 3 of Ref. 41). Because the X-line is located very 
close to one end of the current sheet, the tension force is 
directed predominantly towards the other direction (see 
also Refs. 1231 and 1361). As in previous simulations of recon- 
nection with asymmetric inflo w 23 i 45 i 47 ~ — and asymmet- 
ric outflow ^ 23 ! 29 the X-line and flow stagnation point are 
separated by a short distance. This separation appears 
to be a ubiquitous feature of asymmetric reconnection. 

The X-line retreats at speeds of ~0. 02-0. 06 and the 
flow stagnation point retreats at speeds of ~0. 03-0. 07. 
This is slower than the X-line retreat speed of ~0.1 ob- 
served in Ref. 29. Early in time, the plasma flow at the 
X-line is in the same direction as the rate in change of 
position of the X-line. While the reconnection process is 
still developing, however, the relative positions of the X- 
line and flow stagnation point switch. As a consequence, 
late in time the plasma flow at the X-line is in the op- 
posite direction of X-line retreat. This switch occurs so 
that the flow stagnation point will be located near where 
the magnetic tension and plasma pressure forces cancel. 

To further understand these results, an expression is 
derived for the rate of X-line retreat. In the assumed ge- 
ometry, the X-line retreats in the direction of increasing 
reconnection electric field strength. In the resistive MHD 
limit, the X-line retreats due to either advection by the 
bulk plasma flow or by diffusion of the normal component 
of the magnetic field. 

Interestingly, previous PIC simulations show that 
while the ion flow at the X-point is in the same direction 
as X-line retreat, the electron flow at the X-point is in the 
opposite direction (see Fig. 1 of Ref. |29|) . Because in Hall 
MHD the magnetic field is frozen into the electron fluid 
rather than the bulk plasma, this result suggests that a 
similar effect occurs in fully kinetic simulations. It will 
be important in future work to determine whether or not 
this class of behavior occurs in more realistic geometries 
such as CME current sheets and the Earth's magnetotail, 
or for substantially different plasma parameters. 

The simulation presented in this paper can be used 
to assess the validity of the assumptions made by the 
steady-state model of reconnection with asymmetry in 
the outflow direction presented in Ref. [l|. In particular, 
two assumptions require refinement. The first assump- 
tion is that the magnetic tension force contributes evenly 
to outflow from both sides. In the simulation, however, 
the tension force is almost entirely directed towards the 
unobstructed downstream region because the X-line is 
located very close to the obstructed exit. The second 
assumption is that the current sheet has approximately 
uniform thickness along the outflow direction. However, 
Figs. [Tfc) and [2](a) show that the current sheet has a 
characteristic single wedge shape. Consequently, mod- 
els for asymmetric outflow reconnection will need to be 
refined to account for these effects. 

The results of this paper may have several implications 
for current sheets that form in the wakes of CMEs. In 
these events, the locations of the predominant X-line and 



flow stagnation point are expected to be located near the 
base of the current sheet (see also Refs. [TH, [10, and l37[). 
Consequently, the mass, momentum, and energy fluxes 
are expected to be greater in the antisunward direction. 
This may provide an explanation for some of the observed 
properties of CMEs: that the masses of CMEs increase 
after leaving the Sun^~— and that CMEs are heated 
even after leaving the flare site^ 9 - - — Additionally, cur- 
rent sheets behind CMEs are expected to be thinner when 
observed at low altitudes and thicker when observed at 
high altitudes (e.g., compare the inferred thicknesses in 
Refs. Hand M). However, more observational and the- 
oretical work is necessary before drawing conclusions. A 
complicating factor is that CME current sheets may be 
susceptible to the tearing 63,64 and plasmoid 38 43 insta- 
bilities. 

Effects associated with competition between multi- 
ple competing reconnection sites are likely to have con- 
sequences on the physics of turbulent reconnection. 
The model for turbulent reconnection by Lazarian and 
Vishniac 65 assumes that outflows from each of the small- 
scale reconnection sites in a large-scale current sheet do 
not significantly affect other reconnection sites. The two- 
dimensional simulations reported in this paper cannot 
fully address this problem. Hence in future work we will 
perform three-dimensional simulations with the initial 
perturbations offset from each other in the out-of-plane 
direction. 

The resistive MHD simulations presented in this pa- 
per are not directly applicable to magnetic reconnection 
in the Earth's magnetotail. There, two-fluid and colli- 
sionless effects not included in the resistive MHD frame- 
work are known to be important. Fully kinetic simula- 
tions such as those presented in Ref. 29 show many of the 
same features as the double perturbation simulation pre- 
sented here. Moreover, the shape and structure of the 
out-of-plane quadrupole magnetic field associated with 
two-fluid reconnection can be greatly affected by X-line 
retreat j 29 i 66 Hence, the NIMROD simulations reported in 
this paper will be extended to include two-fluid effects. 
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